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ABSTRACT 

When  a  metal  matrix  composite  is  cooled  down  to  room  temperature 
from  the  fabrication  or  annealing  temperature,  residual  stresses  are  induced  in 
the  composite  due  to  the  mismatch  of  the  thermal  expansion  coefficients 
between  the  matrix  and  fiber.  A  method  can  be  derived  for  calculating  the 
internal  stresses  appearing  in  Metal  Matrix  composites  of  Al  matrix  with  SiC 
particles  due  to  differences  in  thermal  expansion  coefficients.  Special  attention 
is  paid  to  creep  deformation  in  the  matrix  phase.  The  analysis  shows  that 
considerable  internal  stresses  and  creep  deformation  appear  in  the  composites 
when  subjected  to  cooling. 
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I.  INTRODUCTION 

A.  METAL  MATRIX  COMPOSITES(MMCs) 

Discontinuously  reinforced  Metal  Matrix  Composites(MMCs)  represent  a 
group  of  materials  that  combine  the  strength  and  hardness  of  the  reinforcing  phase 
with  the  ductility  and  toughness  of  the  matrix.  Aluminum  alloys(Al)  reinforced 
with  Silicon  Carbide(SiC)  in  particulate,  platelet,  or  whisker  form  and  fabricated  by 
powder  metallurgy  methods  are  receiving  a  great  deal  of  attention  from  researchers 
and  engineers. 

B.  THE  ANALYSIS  OF  MMCs(Al/SICs)  IN  THE  STABLE  MEMBER 
Previous  research  has  evaluated  the  dimensional  stability,  and  the  thermal 

and  mechanical  properties  of  several  Al/SiC  MMCs  in  stable— member  applications 
for  missile  inertial  guidance  system.  The  results  reveal  that,  although  the  candidate 
materials  from  a  powder  blend  of  SiC  and  Al  alloy  consolidated  by  Vacuum  Hot 
Pressing! VHP)  into  cylindrical  billets,  followed  by  Hot  Isostatic  Pressing(HIP)  for 
full  densification.  have  better  isotropic  properties,  these  MMCs  show  some  micro 
creeping  behavior  in  service.  The  microcreep  of  these  MMCs  will  affect  the 
dimensional  stability  of  stable  members,  and  is  likely  to  arise  from  two  sources:  (1) 
creep  of  the  metal  matrix  caused  by  internal  stresses  creep  conditioned  by  externally 
applied  stresses  and/or  (2)  phase  transformations  during  the  creep  condition.  A  first 
step  toward  understanding  the  cause  of  the  dimensional  stability  problem  is  to 
analyze  the  influence  of  internal  stresses.  The  internal  stresses  are  thermal  residual 
stresses  when  a  Metal  Matrix  Composite  is  cooled  down  to  room  temperature  from 
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the  fabrication  or  annealing  temperature.  These  residual  stresses  are  induced  in  the 
composite  due  to  the  mismatch  of  the  thermal  expansion  coefficients  between  the 
aluminum  alloy  matrix  and  the  silicon  carbide  particles.  The  model,  based  on 
Eshelby's  for  mismatch  problems  with  simplified  linear  elastic  material  behavior  for 
both  particles  and  matrix,  has  been  used  to  solve  the  problem  of  thermal  residual 
stress  [ref  1]. 

C.      PURPOSE 

The  purpose  of  this  thesis  is  to  calculate  the  thermal  residual  stress  of  Al/SiC 
using  theoretical  methods  and  then  to  determine  what  influence  this  thermal 
residual  stress  has  on  creep  deformation.  In  order  to  do  this,  the  thermal  residual 
stresses  will  be  determined  by  focusing  on  elastoplastic  matrix  and  elastic  inclusions 
rather  than  on  elastic  matrix  and  elastic  inclusions.  Therefore,  the  case  where  both 
phases(  matrix  and  inclusions  )  are  perfectly  elastic  will  be  treated  first.  Then,  the 
theoretical  results  will  be  compared  with  previously  obtained  experimental  data  for 
Al/SiC  composites. 


II.  LITERATURE  REVIEW 

A.  MMCs 

Metal  matrix  composites(MMCs),  including  eutectic  composites,  are  becoming 
important  in  applications  as  structural  components  which  are  to  be  at  intermediate 
and  high  temperatures.  When  MMCs  are  fabricated  at  high  temperature  or  annealed 
at  certain  high  temperatures,  the  MMCs  have  undesirable  properties,  such  as  low 
tensile  yield  and  ultimate  strengths.  Those  results  are  mainly  due  to  residual 
stresses  that  are  caused  by  the  mismatch  of  the  thermal  expansion  coefficient 
between  the  matrix  and  fiber.  The  residual  stresses  so  induced  have  been  observed 
in  tungsten  fiber/copper  composites  and  in  SiC  whisker/6061  Al  composites  and 
were  based  on  analyzed  a  one— dimensional  model  for  a  continuous  fiber  system  or 
spherical  particle  system.  The  model,  based  on  Eshelby's  equivalent  inclusion 
method,  has  been  used  to  solve  the  problem  of  thermal  residual  stresses[ref  2]. 

B.  ELASTIC  MATRIX  AND  ELASTIC  INCLUSIONS 

Let's  consider  a  two— phase  suspension  of  finite  extent  in  which  one  phase  is  a 

matrix  while  the  other  is  in  the  form  of  inclusions.  It  is  assumed  that  the  inclusions 

can  be  treated  as  identical  ellipsoids  with  corresponding  axes  aligned.  Further,  the 

inclusions  are  assumed  to  be  randomly  distributed  in  the  matrix  in  such  a  way  that 

the  suspension  is  homogeneous  on  a  macroscopic  scale.  The  task  of  the  present 

section    is    to   solve    Eshelby's    transformation    problem    for    the    case    where    the 

fractional  volume  f  of  the  inclusions  is  finite.  Following  Eshelby  [ref  2],  we  consider 

first  an  auxiliary  problem,  what  is  the  resulting  stress  a    in  an  inclusion  when  every 

ij 

inclusion,  having  isotropic  plastic  moduli  equal  to  those  of  the  matrix,  is  subjected 


T 
to  a  uniform  transformation  strain  €  To  find  the  exact  solution  appears  to  be 

i] 
impossible  since  no  detailed  information  is  given  about  the  spatial  distribution  of 

the  inclusions.  All  we  may  find  in  this  situation  are  approximate  average  fields.  The 

average  constrained  strain  €  ~  in  anv  inclusion  may  be  written  in  the  form 

ij 

i  n 

6"°  =  6  C  +  6  C  (  1  ) 

i 

c 

Here,  €        is  the  constrained  strain  in  a  typical  inclusion  which  would  be  produced  if 
ij 

it  alone  is  embedded  in  an  infinite  elastic  body.  This  term  is  given,  according  to 
Eshelby[ref    ],  as 


€  C    =  S        €  T  (  2  ) 

ij  ijkl     ij 


Where  S        is  Eshelby's  tensor  which  depends  on  the  aspect  ratios  of  the  ellipsoidal 
ijkl 

inclusion  and  Poisson's  ratio. (Here  and  in  the  following  we  use  the  usual  summation 

ii 

c 

convention:  the  range  of  subscripts  is  form  1  to  3).  6"      is  a  contribution  due  to  all 

i  J 

the  surrounding  inclusions  and  the  presence  of  free  boundary;  this  term  is  entirely 

elastic,    and    is   assumed    to   be   eiffectively    constant    in   and    around    the    typical 

1 1 

c 

inclusion.  Then,  it  will  be  a  good  approximation  to  regard  €"  as  the  average 
elastic  strain  in  the  matrix  for  the  contribution  of  any  one  inclusion  to  the  average 
field;  outside  it  is  negligible  when  there  are  a  great  number  of  inclusion  in  the 
material.  For  a  homogeneous  elastic  body  in  static  equilibrium,  a  volume  integral  of 
any  component  of  the  elastic  strain  associated  with  an  internal  stress  should  vanish 
provided  that  the  integration  is  taken  throughout  the  body.  Thus 


II 

(i— f)  €"c  +  f(e~c-e~T)  =  0  (  3 


From  (1),  (2)  and  (3),  we  have 


€"°  =  (l— f)S     eT  +  feT  (4) 

i  j  ljkl    k  1  i  j 


The  average  internal  stress  in  any  inclusion  is  given  by 


al  =  C      (€-G-€T)  (5) 

ij  ijkl       kl  kl 


Where  C        is  the  elastic  moduli  of  the  matrix  (and  also  of  the  inclusions  in  the 
ijkl 

present  situation)  [ref  3].  Where  the  inclusions  have  elastic  moduli  differing  from 

those  of  the  matrix,  the  same  approach  as  used  by  Eshelby  [ref  2]  will  be  applicable. 

* 
The  essence  of  this  approach  is  terms  of  equivalent  inclusions.   Let   C         be  the 

ij  kl 

T* 
elastic  moduli  of  the  inclusions,  and  6       be  a  uniform  transformation  strain  of  everv 

1J.  ...  T 

inclusion.  Then  the  transformation  strain  of  the  equivalent  inclusions,  6     ,  can  be 

i  j 

determined  by  solving  a  set  of  equations 


C      (  €-c-e  T)  =  C*    (  6-C-6T*  )  (6) 

ijklv     kl  kl    '  ijkl       kl  kl  '  v 


with  (  4  ).  The  average  stress  in  any  inclusion  is  given  by  (  5  ),  of  course.    Consider 
a  composite  specimen  which  is  free  of  stress  at  a  references  temperature  To-  Let  the 

matrix  and  the  inclusions  have  isotropic  linear  thermal  expansion  coefficients  a  and 

* 
a  ,   respectively.   The  elastic  state  when  the  specimen  is  subjected  to  a  uniform 


temperature  change,  aTo=T— To,     is  described  by  the  situation  where  the  matrix 

and      the      inclusions      undergo     uniform     transformation      strains      qaT8       and 

*  13 

a  aT<5  respectively.  Here  6     is  the  Kronecker  delta.  The  resulting  state  of  stress 

can  be  analyzed  by  putting 


€T*  =  (a*-a)  AT  6  (7) 

** 

The  overall  thermal  expansion  coefficients  a      can  be  expressed  as  a  function  of  a, 

*  * 

a  ,  C       ,  C       ,  f  and  the  aspect  ratios  of  the  inclusions, 
ijkl       ij  kl  K 


C.  ELASTOPLASTIC  MATRIX  AND  ELASTIC  INCLUSIONS 

When  the  temperature  change  is  great,  one  should  expect  plastic  deformation. 

either  in   the  matrix  and/or  in   the  inclusions,   which  relaxes  the  internal  stress 

accumulated.  An  important  practical  case,  the  case  where  the  plastic  deformation 

occurs   only    in    the    matrix,    will    be    considered    in    the    following.    Further,    for 

simplicity,  we  assume  that  the  plastic  strain  is  uniform  throughout  the  matrix  and 

the  plastically  nondeforming  inclusions  can  be    treated  as  identical  spheroids.  Then, 

the  plastic  strain  in  the  matrix,  €    ,  may  be  written  in  consideration  of  a  symmetric 

ij 

mode  of  plastic  deformation  and  the  volume  constant  law,  as: 

6P  =  €  \6   6    -k6  6+6  6))  (8) 

ij  PL   3i   3j        ^     li    lj  2i  2j  J 

** 

Where  6n  is  a  temperature— dependent  parameter.  In  this  case,  the  overall  strain  € 

™  *i 

—  T 

of  the  composite  is  given,  with  use  of  Equation  €    =f€     ,  by 

ij        ij 


**  rp 


€       =  oaTS   +  6  P+f€  i  (  9  ) 

T 
Where  6      is  calculated  from  (4)  and  (6)  with 

£T*  =  (a*  -a)*T6  -€p  (  10  ) 

Now,  the  temperature  dependence  of  6p  will  be  found  through  an  energy— balance 
consideration.  Since  energy  dissipative  process  is  involved  in  plastic  deformation, 
the  decrease  in  the  elastic  energy.  —  6Eei,  corresponding  to  a  virtual  plastic 
deformation  8£  ,  should  be  balanced  against  the  dissipation  energy.  If  the  matrix 
material  has  a  constant  flow  stress  a     in  a  uniaxial  deformation,  the  condition  is 

y 

written  as: 


-ffiei  =  (l-iVy|<%|    Uy>0)  (11) 


Where  Eei  =  —  Acf  € 

-     ij    ij 

=  C,€2  +  C2€p(a   -a)AT 

+  C3(Q*-a)2AT2  (12) 


It  is  noted  that  when  the  composite  is  elastoplatic  the  following  condition  should  be 
satisfied: 


(EeOj^EeOn 


Where  the  subscripts  I  and  II  refer  to  the  entirely  elastic  case  and  the  elastoplastic 
case,  respect ively[ref  3]. 

D.   CREEP 

If  a  stress  is  suddenly  applied  to  pure  metals,  some  solid  solutions  and  most 
engineering  alloys  at  a  temperature  near  or  greater  than  0.5Tmf where  Tm  is  the 
absolute  melting  point  of  the  metal  or  alloy)  deformation  proceeds  as  shown  in  Fig. 
2.1.  The  initial  application  of  stress  causes  an  instantaneous  elastic  strain  €o  to 
occur.  If  the  stress  is  sufficiently  high  art  initial  plastic  deformation  6P  also  occurs. 
At  low  temperatures,  significant  deformation  ceases  after  the  initial  application  of 
stress  and  an  increase  in  stress  is  required  to  cause  further  deformation.  At  elevated 
temperatures,  deformation  under  a  constant  applied  load  continues  with  time.  The 
early  stage  of  such  deformation  called  primary  creep  is  characterized  by  an  initially 
high  creep  rate  d€/dt  which  gradually  decreases  with  time.  Eventually  a  linear 
variation  of  creep  strain  accumulation  with  time  is  observed.  This  steady— state 
creep  region  is  characterized  by  a  constant  minimum  creep  rate.  Steady— state  creep 
rates  depend  significantly  on  stress  and  temperature  and  are  used  frequently  to 
compare  the  creep  resistance  among  alloys.  After  significant  deformation  in 
steady— state  creep,  necking  occurs  or  sufficient  internal  damage  in  the  form  of  voids 
or  cavities  accumulate  to  reduce  the  cross— sectional  area  resulting  in  an  increase  in 
stress  and  creep  rate.  The  process  accelerates  rapidly  and  failure  occurs.  This  region 
of  the  creep  curve  is  called  tertiary  creep.  Fig.  2.1  also  shows  the  derivative  of  the 
creep  curve  or  the  creep  rate  curve.  Fig.  2.1  represents  the  most  common  creep 
strain/time  behavior  and  Fig.  2.2  illustrates  the  strain  rate/time  behavior  at 
temperatures    below    approximately    0.3Tm.    In    Fig.    2,    only    transient    creep    is 


observed    which    is   characterized    by    a    continuously    decreasing   creep    ratt-    that 

approaches  zero  as  the  inverse  of  time.   Such  low— temperature  creep  behavior  is 

called  logarithmic  creep.  Fig.  2.2  also  illustrates  the  creep  rate/time  behavior  of 

alloys  that  exhibit  a  continuously  increasing  creep  rate  in  the  early  stages  of  creep. 

Such  alloys  have  been  designated  Class  1  alloys.  Creep  by  viscous  dislocation  glide 

results  when  dislocations  glide  or  move  in  slip  planes  under  the  action  of  an  applied 

stress.  These  dislocations  drag  along  solute  atoms  attracted  to  the  strain  fields  of 

the  dislocations.  In  order  for  the  dislocations  to  move  the  solute  atmospheres  must 

diffuse  in  the  direction  of  dislocation  motion.  At  lower  temperatures(-0.3Tm)  or  at 

higher  stress  levels,  creep  occurs  by   thermally  activated  dislocation  glide.  Under 

these  conditions  dislocations  can  overcome  barriers  to  motion  without  dislocation 

climb.  The  localized  motion  of  short  segments  of  small  dislocations  is  important  in 

overcoming  barriers.    Because   of  the   larger  contribution   of  stress   in    thermally 

activated  dislocation  glide  a  temperature  and  stress  dependence  for  creep  different 

from  that  for  dislocation  climb  or  viscous  glide  will  be  observed. fref  6] 
Cneep  rate,  de/dt  Creep  strain.  £ 
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Fig.  2.1  Variation  of  creep  and  creep  rate  with  time 
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Fig.   2.2  Creep  rate  vs  time 
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Fig.   2.3  Typical   creep  curve 
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III.  THEORETICAL  BACKGROUND  AND  DERIVATION 

A.   THEORETICAL  BACKGROUND 

To  simulate  Metal  Matrix  Composites(MMCs),  we  consider  an  infinite 
isotropic  body  containing  an  ellipsoidal  fiber  as  shown  in  Fig.  3.1.  The  infinite 
body( composite)  experiences  a  temperature  Tl  —  TO.  where  Tl  is  assumed  to  be 
larger  than  the  relaxation  temperature  of  the  matrix  Tr  and  of  course,  Tr>T0.  To 
make  the  problem  simpler,  the  case  of  one  fiber  embedded  in  the  infinite  body  is 
considered  here. 


Fig.  3.1.  Theoretical  model 
Based    on    the    model    in    Fig    2,    let     us    consider    the    case    of    elastic 
inclusions— elastic     matrix,     elastoplastic     matrix— elastic     inclusions     and     creep 
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deformation: 


1.        Elastic  Inclusions  and  Matrix 

-C 
First   let  us  consider  the  average  internal  constrain  strain  €     .  To  find 

-C  .  ij 

6..,     we     insert     Equation(l)     and     Equation(3)     into     Equation(4)     and     find 

-C  T  T         • 

€     =(1— f)S       €       +  f€     ,  which  represents  the  average  internal  constrain  in  anv 

i  j  ij  kl    k  1  i  j 

inclusion  in  terms  of  transformation  strain.  Next,  the  average  internal  stress  in  any 

T 
inclusion   can   be   represented    by   Equation(5).    Here,    €       is   an    unknown   term. 

ij 

Therefore,  to  find  this  unknown  term  we  must  consider  Equation  (6)  and  Equation 

T* 
(4).   In  Equation  (6),   *   represents  the  inclusion,  and  G        represents  the  uniform 

1]         .  T 

transformation  strain  for  every  inclusion.  Once  we  have  determined  €     ,  according 

T  1J    • 

to  the  Equation  €    =  f€       we  can  determine  the  overall  strain  of  the  specimen,  and 

ij  i  j  m. 

T* 
in  accordance  with  Equation(7)  we  can  determine  the  uniform  transformation  6 

>  j 


2.        Elastoplastic  Matrix  and  Elastic  Inclusions 

The  case  for  elastoplastic  matrix -elastic  inclusions  is  different  from  that 

of  elastic  inclusion  and  matrix.  This  is  due  to  the  fact   that  in  this  case  plastic 

deformation  occurs  only  inside  the  matrix,  so  we  must  consider  the  plastic  strain  in 

the    matrix.    Therefore    we    can    find    €  ^    in    accordance    with    Equation(S)    and 

ij 

Equation(lO).     From    Eei     =    —    Av   £         we    can     represent     a'  =F(€  *j    and 

*    ij    ij  ij  ij 

T  D 

€      =F(€  p).  After  representing  them  this  way, 
i  j  ij 

Eel  =  -  A  O    6 

*J         *    J  r.  ~ 

2     2 

=C1  €  2  +  C2  €  (a*  -  a)AT  +  C2  6  (a*  -  af  aT  (  14  ). 

p  P  P 
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Differentiating  this  with   respect   to  G   ,  we  can  obtain  the  value  of  G     as  shown 
below: 


C2      *  1 

6p  =  ~2UT  ( a    _  °)aT  Z  7TTT^ l  ~fK 


+   2( 


(15) 


Therefore,   using  the  value  of  6     to  find  the  value  G  ^,   the  value  of  a'    can  be 

P  ij  ij 

obtained,    and    inserting   this   value   into   Equation(16),   we   are   ultimately   able   to 
obtain  the  total  average  stress  in  the  matrix. [ref  9] 


<  *A.  >™  =  C^  <  ekl  >Vr  =  -f  J 


i] 


i) 


(16 


3.        Creep 

Applying  the  value  of  the  internal  stress(above)  to  Equation(13).  we  can 
determine  creep  deformation.  Following  from  this  theoretical  base,  let  us  first 
calculate  mathematically  the  thermal  residual  stress  of  Al/SiC  MMCs. 


B.      DERIVATION 

1.        Elastic  Inclusions  and  Elastic  Matrix 

T 
First  let  us  calculate  G      for  the  case  of  elastic  inclusions  and  elastic  matrix. 

T      U 
In  order  to  determine  G     ,  we  rewrite  Equation(4)  in  matrix  form: 

1 J 


-c 

€!c 

Eg 

G66 


=(l-f) 


S  1 1  S  12  S  1 3  S  14  S  15  S  1  6 

€i' 

61 

S21  S22  S23  S24  S25  S26 

62: 

€2 

S31  S32  S33  S34  S35  S36 

63! 

+  f 

€3 

S41  S42  S43  S44  S45  S46 

e4! 

€4 

S51  S52  S53  S54  S55  S56 

e5! 

e5 

.  Sei  Se2  Se3  Se4  Ses  Se6 . 

Eel 

.U 
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Simplifying  this  matrix  we  obtain: 


C2') 
63£ 

C66 


=(l"f) 


€lt 


Suf       +     S,2 

S2l€|j    +€2$(    S?2 


S  3 1  e  ii  +S32  62$  +  €s§( 
644(244   + 
65KS55   + 
€ee(  See    + 


6i|  +    S,3€33 
+  i   )    +  S  2  3633 
S33  +  f  ) 


f 

n 


r6,u 

€0^ 

+  f 

€44 

Gsj 

.€6fc. 

Further,  simplifying  this  matrix  we  find: 


€if  =  €,Hl-f)S„  +  €2|SI2(l-f)  +  esiS,3(l-f)  +  fe,t 
€^2  =  €,Ki-f)S2.  +  €2^22(1-0  +  63^23(1-0  +  teA 
€33  =  6.K1-OS3.  +  €2^32(1-0  +  €3^33(1-0  +  fe3& 
€iS  =  64i(l~f)S44  +  iu\ 
eis  =  6bK1-I)Sbs  +  fe5| 

666  =  666(l-f)S66  +  ^66 


(17) 


Refer  to  Appendix  (A)  for  the  values  of  Sijkl- 
Next,  let  us  rewrite  Equation(6)  in  matrix  form. 


G 


H 


(l-i/)    v       v       0     0     0' 

* 

G2g 

€33 

€ij 

" 

v    (l-i/)    i/000 

€2| 

v       v    (l-i/)    0     0     0 

€3! 

0       0       0     2/i    0     0 

^ 

€4$ 

0       0       0       0  2/i    0 

ej£ 

€5$ 

0       0       0       0     0  2/i  _ 

- 

Egg 

.€68. 

- 

(l-i/)    1/       1/      0     0     0' 

- 

[€iil 

€1  * 

1/    (l-i/)    1/      0     0     0 

€22. 

€2J* 

1/       1/    (l-i/)    0     0     0 

€33 

€3S* 

0       0       0     2/i    0     0 

€<* 

€4}* 

0       0       0      0  2/j    0 

€5£ 

€sf* 

0       0       0      0    0  2/i 

€66  . 

€6| 
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Simplifying  the  above  matrix,  we  have  the  following: 


i.*  ►* 

-i/He25  -^He33 

(l-i/)H62^*-  */He3f 

eni/J+€«i/J+€§3(l-i/)J=i>G€if+  i^€2|+(l-i/)G€3^-i/H6,$  - 

i/H62|*-(l-^H€3f 

--C 


* 


eee  =  0  (18) 


T 
From  Equation  (17)  and  (18),  if  we  simplify  and  solve  for  €     ,  we  obtain: 

ij 


*  *  * 

Keii+Le2$+M€3S=-(i-^Heii  -vtt^  -j/H€3& 
Ne,HP€2HQ€3^=-t/He1f-(l-i/)He2^-*/H€3S* 
R€,i+Se2l+T63h-fHe,i*-f/He2^-(l-»/)He35* 
€44  =  €44 

€«J  =  £*C  (  19  ) 


Now,  we  have  6  unknown  values  and  6  Equations.  Therefore,  applying  Equation(19) 

T 
to  determine  € 


.     -YL-XM  -HA(l+>) 
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e4S  =  e4S 
e6£  =  €6fc 


W  -  vx 
— 0 

ZU-XW 
UX  -  XV 


(20) 


Putting  the  above  equations  into  simplified  form: 


[€ijl 

r  z,  1 

623 

Y, 

€T  = 

e3& 

€44 

= 

C4  4* 

e5J 

€ftU 

[€68  J 

ltd  J 

Here,  refer  to  Appendix(C)  for  constants  B  through  Z.  Thus,  inserting  the  value 

T 
€       into  Equation(5),  which  is  the  average  internal  stress  in  any  inclusion,  we  can 

M 
determine  the  internal  stress. 


.  I 

^33j 
^4  4  j 

<766 


=  G 


'(1-1/)      V          V          0  0  0 

v    (l-i/)    v       0  0  0 

v    '  1/    (l-i/)    0  0  0 

0       0       0     2/*  0  0 

0       0       0       0  2p  0 

0       0       0       0  0  2/t 


BZ'+CY'+DX'-Z1 
CZ'+BY'+DX'-Y' 
FZ'+FY'+EX'-X1 

0 

0 

0 


2.        Elastoplastic  Matrix  and  Elastic  Inclusions 

Up  to  this  point  we  have  found  the  internal  stress  in  the  case  of  an 
elastic  matrix  and  elastic  inclusions.  Now  we  will  consider  the  case  of  elastoplastic 
matrix  and  elastic  inclusions.  Here  also  the  method  of  determining  the  internal 


16 


and  elastic  inclusions.  Here  also  the  method  of  determining  the  internal  stress  is 
similar  to  the  case  of  elastic  matrix  and  inclusion.  However  in  this  case,  since  6    , 

the  plastic  strain,  is  included  in  the  term  6       ,  the  next  case  is  more  complicated. 

.lj  d  T* 

Therefore,  before  we  consider  the  relationship  between  €  v  and  €       ,  let  us  first  look 

T* 
at  €       : 

i  j 

€  T*=(a*  -  a)AT6  -€  [8   8    -1(8  8  +  6  8  )]  (  6  :  i=j-»l,  tfH  )•  Rewriting  the 

ij  ij       P     3i   3j        *•       li    lj         2i  2j  ij 

above  equation  in  matrix  form  we  have 


e3 
€4 

€51 
€6 


=(a  -a)AT 


■*li' 

€"d 

^22 

e22P 

^33 

€33? 

^44 

6«£ 

^55 

€56P 

.  ^66  . 

€66P 

(21) 


But  according  to  Equation(8),  €  ^  can  be  shown  as  follows: 

ij 

£31  £31— 1/2(  8\\  £11  +  621  $21 

^12  ^12— 1/2(  ^12  ^12+^2  ^22 

^33  ^33— 1/2(  $33  ^33+^23  ^23 

634  634— l/2f  £|4  #14+ #24  ^24 

^35  ^35— 1/ A  ^15  ^15+^25  ^25 

^36  ^36— I/A  ^16  ^1 6+^26  ^26 


re"ni 

e< 

=€ 
P 

.  666P  . 

From  the  above  equation  we  can  obtain  the  following: 


-1/2€P  1 

T\ 

-1/2€P 

€P  = 

eP 

»j 

0 

0 

0 

(22) 


Therefore,  inserting  £quation(22)  into  Equation(2I)  and  simplifying  we  obtain: 
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6  J*=(a*  -  a)AT  +  €p/2 
eJ*=(a*-a)AT  +  €p/2 
«„   =(<■   -«>T-6p 

eT*=0 

44 

€T'=0 

55 

€  T*=0  (23) 

66 

Next,  using  Equation(4)  and  Equation(6),  if  €      is  expressed  as  ETi,  the  result  is  the 

\j  T 

same  as  Equation(17).  Inserting  Equation(17)  into  Equation(6)  and  expressing  € 

T*  .  ij 

as  £        we  obtain: 

rp  rp*  rp*  rp* 

€      =Qi€i  |     —  Q2^22    —  Q3^33 
T  T*  T*  T* 

e    =P4en  +P5622  +P6€ss 
2  2 

T»  T*  T*  T* 

6   t=Pi€ii    +P2€22    +P3€33 
33 

6T=0 

44 

ej=0 

Ej-0  (24) 


c  T  • 

Next,  in  order  to  find  €11,  we  insert  the  value  of  €      into  Equation(14)  and  simplify: 

»i 

_c  »j»*  »p*  'p* 

€|l=R|€||    +  R2€2  2    +  Rs€3  3 


€22=FU€|  1    +  R5G22    +  R6€3  3 
€SS=R7€|  I    +  R^2  2    +  R«€3  3 


€i4=0 
€g5=0 
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€66=0 


(25) 


Therefore  a\\  can  be  found  in  accordance  with  Equation(5).  (The  value  of  this  o~^  is 
different  from  the  value  of  <rij  for  elastic  inclusion  and  elastic  matrix,  the  reason 
being  that  <?{)  for  elastic  inclusion  and  elastic  matrix  does  not  contain  €..).  Here,  in 


-I 


ij 


order  to  determine  the  value  of  a^  we  must  first  find  the  value  of  €  .  From 
Equation(lO),  if  we  differentiate  Eei  with  respect  to  €  we  can  obtain  Equation(ll). 
Therefore,  from  ffid  =[  2C|€  +  C2(a  -  a)AT]fe  and  -6E4  =  (l-i)<r  \fe  |,  we 
can  obtain  the  value  of  €  .  Again,  inserting  this  value  of  €  into  Equation(22)  we 
can  find  the  value  of  £  *\  Finally,  by  inserting  the  value  of  G  ^  into  equation(27)  we 
obtain  the  value  of  aij.  Utilizing  Equation(lO)  again  to  find  the  value  of  € 


Ed=-if 


-  In 

^66 


* 

e2f* 

£3§* 

€5|* 
€66 


-¥ 


X,A'  +  X2€p 
X3A1  +  X4€p 
X5A'  +  X6GP 

0 

0 

0 


A'  +  €p/2 
A'  +  €p/2 
A' -€p 

0 

0 

0 


X? .  Xi 


=  -  h  (eP(  f+¥  -  X«  )  +  €pA'(  7i+XJ+^-X4-X8+X«  ) 
+  A,2(  X,+X3+X5 )] 


(26) 


Therefore,  differentiating  Eei  in  respect  to  €p: 


ffiel  =-t(  C2A'p+  26pC,  )& 


(27) 


And  inserting  Equation(27)  into  Equation(ll)  and  simplifying  obtains: 


f(  C2A'p+  2GpC,  )&=(l-f)ay 


(28) 
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If  we  find  the  value  of  €p  from  Equation(28). 


€p=-  §^( a   -  a)AT  +  i^(l-f)ay  (  27  ) 

(  +  :  fe  <0  -»  Cooling,  -:  fc  >0  -  Heating  ) 


T 
If  this  value  of  €p  is  inserted  into  Equation(24)  and  (25),  the  value  of  €      and  the 

value  of  €     can  be  found,  and  if  these  values  are  inserted  into  Equation(5)  the  value 

-I         1J 
of  a     can  be  determined.  Therefore,  inserting  this  value  into  Equation(16)  we  are 

ultimately  able  to  obtain  total  average  stress  in  the  matrix.  Refer  to  Appendix(  C  ) 

for  the  constants  and  actual  values. 


3.        Creep 

Up  to  this  point  we  have  been  determining  the  internal  stress  both  in  the 
case  of  elastic  matrix  and  elastic  inclusion  and  in  the  case  of  elastoplastic  matrix 
and  elastic  inclusion.  But  in  these  two  cases  the  actual  effect  on  creep  deformation  is 
the  internal  stress  in  the  case  of  elastoplastic  matrix  and  elastic  inclusion. 
Accordingly,  by  inserting  the  value  of  the  internal  stress  into  Equation  (13),  we  can 
plot  the  microcreep  deformation  phenomenon  as  a  function  of  time  as  shown  in  Fig. 
3.2.  Refer  to  Chapter  IV  part  A  for  more  detailed  information. 
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FIG.  3.2   MicroCreep  as  a  function  of  time  t 
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IV.  RESULTS  AND  CONCLUSIONS 

A.      RESULTS 

The   thermo— mechanical   data  of  the   matrix   and   fiber   for   the   theoretical 
calculations  are  obrtained  from  the  [ref  13]. 
Annealed  2024  Al  matrix: 
Em  =  47.5  GPa 
a    =  47.5  MPa 

y 

v  =  0.3 

a  =23.6xlO"6/K  (31) 

SiC  fiber: 

Ei  =  427  GPa 

f  =  0.2 

a*  =  4.3X10-6/  K 

l/d=  1.5  (32) 

Where  the  average  value  of  the  fiber  aspect  ratio(l/d)  was  used  [ref  13].  The 
temperature  drop  aT  is  define  as 

aT  =  T1-T0  (33) 

Where  Tl  is  taken  as  the  temperature  below  which  dislocation  generation  is 
minimal  during  the  cooling  process  and  To  is  the  room  temperature.  Thus,  for  the 
present  composite  system  aT  is  set  equal  to  — 200K.  From  the  data  given  by 
Equation  (31),  (32)  and  the  use  of  <  a  >m  =  — ia  ,  we  have  computed  the  stresses. 
Next,  the  thermal  residual  stresses,  averaged  in  the  matrix  of  SiC  fiber/2024  Al,  are 
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predicted  by  (17)  and  the  result  on  <  a~    >a  as  plotted  in  Fig.6,  where  33  denote 

*  Jo 

the  component  along  the  longitudinal  direction(z).  The  average  theoretical  thermal 
residual  stress  is  predicted  to  be  tensile  in  nature,  and  the  average  residual  stress  in 
the  longitudinal  direction  to  be  larger  than  the  average  residual  stress  in  the 
transverse  direction  [Table  5].  The  fiber  aspect  ratio(o=l/d)  of  SiC  fibers  has  been 
observed  to  be  variable  [Appendix  A,B]-  In  the  present  model  we  have  used  the 
value  of  1/d,  1.5  to  predict  the  thermal  residual  stress  of  the  composite. 


TABLE  1:  The  value  of  6  p  (a=l,1.5) 


a(l/d) 


1 

1.5 


6  P 
l  i 


-0.0017 
0.0014 


-0.0017 

0.0014 


€  p 

33 


-0.0017 
-0.0026 


TABLE  2:  The  value  of  6*  j  (a=l,  1.5) 


a(l/d) 


1 
1.5 


6T* 
i  i 


0.0056 
0.0025 


7F*" 


0.0056 
0.  0025 


€T* 

3  3 


0.  0056 
0.  0066 
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TABLE  3:  The  value  of  e~j  |  <v=l,  1.5; 


Q(l/d) 

6"^ 

i  i 

0  «> 

33 

1 
1.5 

0.0053 
0.0025 

0.0053 
0.0025 

0.0053 
0.0060 

TABLE  4:  The  value  of  €;}  (o=l,  1.5) 


3  3 


cKl/d 


1 

1.5 


i  i 


0.0076 
0.0025 


0.0076 
0.0025 


0.0076 
0.0115 


TABLE  5:  The  value  of  a'1  (o=l,  1.5) 


ofl/d) 

<r](MPa) 
1 1 

^-i(MPa) 

22 

^(MPa) 

33 

1 

1.5 

-276.2  3 
-154.82 

-276.  2  3 
-154.8  2 

-276.  2  3 
-154.8  2 

TABLE  6:  The  value  of  «r~l>m  (o=l,  1.5) 


a(l/d) 


1 
1.5 


<<r-i>(MPa) 
1 1 


55.245 
30.905 


«7~i>(MPa) 

22 


55.24  5 
30.905 


<fr~1>(MPa) 

33 


55.245 
70.998 
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TABLE  7:  Comparison  value  of  <  tr     >m  with  TODAY'S 


o(l/d) 

OUR  MODELS  '  S 

TAYA ' S 

<  dj\    >D(MPa) 

<  vl\   >n(MPa) 

1.8 

71  .  156 

67.  894 

As  seen  in  tables  1.  2,  3,  4,  5,  6,  and  7  we  can  observe  that  as  a  increases,  the  values 


P     t* 


-c 


T     -I 


-I 


of  6ij ,  €  i  j,  €ij,  6i  j ,  ^Tj,  and  <  <t7j  >m  also  increase 
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Also  in  Fig.  4.1  we  observe  that  when  the  aspect  ratio(a)  increases,  the  value  of  the 
internal  stress  also  increases  with  it. 


Fig.  4.1   Internal    stress  vs  aspect  ratio 
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And  in  Fig.  4.2  we  can  also  see  that  the  value  of  the  internal  stress  increases  with  an 
increase  in  temperature. 


Fig.  4.2   Internal    stress  vs  temperature  change(&T) 
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In   Fig.   4.3,   in   the  case  of  an  ellipsoid  (a=1.5),   the  values  of  internal  stresses 
increase  with  an  increase  in  volume  fraction, 


Fig.   4.3   Internal    stress  vs  volume   fraction(-f) 
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Next,  if  we  compare  the  theoretically  predicted  value  of  the  internal  stress  with  the 
value  determined  by  R.  J.  ARSENAULT  and  M.  TAYA,  we  see  that  our  model's 
value  is  larger  than  the  value  which  is  obtained  by  R.   J.  ARSENAULT  and  M. 
TAYA,  as  shown  in  Table  7.  They  used  the  material  properties  as  follows: 
Annealed  6061  Al  matrix: 

E=47.5  GPa 

a  =47.5  MPa 

y 

v  =0.33 

a  =23.6xlO"6/K 
SiC  whisker: 

E1=427GPa 

f  =0.17 

q*=4.3x10^/K 

l/d=l.S 
Finally,  Fig.  3.2  is  a  graph  showing  the  strain  as  a  function  of  time,  when  the  aspect 
ratio    is    1.5    and    the    value    of    the    thermal    residual    stress    is    inserted    into 
Equation(13).    If  we   analyze    Fig.    2.3    we   can   see    that    the   creep   deformation 
phenomenon  is  due  to  the  internal  thermal  residual  stress  of  the  MMCs. 

B.      CONCLUSIONS  AND  RECOMMENDATIONS 

The  object  of  this  research  is  to  obtain  the  value  of  the  thermal  residual  stress 
of  an  Al/SiC  composite  using  Eshelby's  theoretical  model,  and  then  to  determine 
what  effect  this  thermal  residual  stress  value  produces  on  creep  deformation. 
Because  exact  creeping  behavior  of  this  material  is  difficult  to  determine  when 
analyzing  creep  deformation,   Andrade's  model   was   used   with   properly   selected 
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constant  values.  But  through  this  theoretical  approach  the  following  conclusions  can 
be  obtained. 


1.  Thermal  residual  stress  due  to  the  difference  in  thermal  expansion 

coefficients  may  be  estimated. 

a.  By  means  of  the  theoretical  model,  we  can  see  that  the  thermal 

residual  stresses  increase  when  the  volume  fraction  of  the  inclusions 
increases. 

b.  In  the  case  of  sphere  inclusions,  that  is,  when  the  aspect  ratio  q=1,  the 

lateral  stress  is  equal  to  the  longitudinal  stress. 

c.  In  the  case  of  an  ellipsoid  inclusions  where  a=l.b  the  longitudinal  stress 

is  greater  than  the  lateral  stress. 

d.  The  thermal  residual  stresses  of  Al/SiC  composites  increased  when 

the  value  of  the  aspect  ratio  of  the  inclusions  increases. 

2.  Microcreep  deformation  can  be  estimated  in  the  model  by  using  the 

thermal  residual  stresses  and  the  Andrade's  model  of  creep 
deformation 

a.  Microcreep  deformation  is  due  to  the  internal  thermal  residual 

stress  of  the  MMCs. 

b.  Dimensional  stability  of  components  will  be  influenced  by  the 

behavior  of  the  composite. 

3  Recommendations 

a.  Presently,  this  thesis  has  only  dealt  with  average  thermal  residual 

stress  from  a  overall  point  of  view,  and  we  need  more  detailed 
local  residual  stresses  surrounding  SiC  particles  should  be 
estimated  for  analyzing  creeping  behavior. 

b.  In  the  future,  research  should  be  done  concerning  relaxation  during 

cooling. 

c.  The  creep  behavior  used  in  the  current  model  is  the  Andrade's 

approximation.  The  real  creep  deformation  for  this  Al/SiC 
composite  should  be  further  studied  in  order  to  get  better 
understanding  of  the  microcreep  of  the  material  based  on 
approximation  values. 
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APPENDIX  A 
PROGRAM  FOR  VALUE  OF  S^o^LS  -  5) 
1.  PROGRAM(FORTAN) 


REAL* 8  NU, ALP, G, TMP, Sllll, S3333,S1122,S1133,S3311 
NU-0.3 

WRITE (6 , 90)  i 

WRITE (6, 95) 
ALP-1.5 
300   TMP-ALP**2-1. 

G- (ALP* (ALP*DSQRT (TMP) -DLOC (ALP+DSQRT (TMP) ) ) ) / (TMP** (3 . /2 . ) ) 
S1111«(3.*ALP**2)/(8.*(1.-1IU)*TMP) 
++(l.-2.*NU-9./(4.*TMP))*G/(4.Ml-NU)) 

S3333-((1.-2.*NU+(3.*ALP**2-1.)/TMP-(1.-2.*NU+3.*ALP**2/TMP) 
+*G))/(2.*(1.-NU))  T  r 

S1122«(ALP**2/(2.*TMP)-(1.-2*NU)-3*G/(4.*TMP))/(4.*(1.-NU)) 
S1133— ALP**2/ (2 . * (1 . -NU) *TMP) + (3 . *ALP**2/TMP- (1. -2 . *NU) ) *G 
+/(4.*(1.-NU))  i 

S3311— (l.-2.*NU+l./(TMP))/(2.*(l'.-NU))  +  (l.-2.*NU+3./(2.*TMP)) 
+*G/(2.*(1.-NU))  ' 
WRITE ( 6 , 100 )  ALP , Sllll , S3333 , S1122 , S1133 , S3311 
ALP-ALP+.l  I 

IF (ALP. GT. 5.)  GOTO  400 
GOTO  300 
90    FORMAT(lX,'      ALP        Sllll       S3333 

+  S3311') 
9  5    FORMAT  ( IX ,  •    ——■—————■■——. 
♦««« -•) 

100   FORMAT (IX, 6 (F10. 5)) 
400   STOP 
END 


S1122     S1133 
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2.  S..v,  VALUE  WITH  ASPECT  RATIO 
ijkl 


ALP 


sim 


S3333 


S1122 


S1133. 


S3311 


1.50000 

.58078 

.39378 

.01421 

.08396 

.01945 

1.60000 

.58828 

.37371 

•01603 

.08992 

.01623 

1.70000 

.59495 

.35518 

.01764 

.09549 

.01350 

1.80000 

.60089 

.33804 

.01906 

.10071 

.01117 

1.90000 

.60620 

.32215 

.02032 

.10560 

.00918 

2.00000 

.61098 

.30739 

.02145 

.11018 

.00747 

2.10000 

.61529 

.29366 

.02246 

.11448 

.00600 

2.20000 

.61920 

.28087 

.02337 

.11851 

.00472 

2.30000 

.62274 

.26893 

.02419 

.12231 

.00362 

2.40000 

.62598 

.25777 

.02493 

.12588 

.00266 

2.50000 

.62893 

.24732 

.02560 

.12924 

.00183 

2.60000 

.63164 

.23752 

.02622 

.13241 

.00111 

2.70000 

.63412 

.22832 

.02678 

.13540 

.00048 

2.80000 

.63641 

.21967 

.02729 

.13822 

-.00007 

2.90000 

.63852 

.21152 

.02776 

.14090 

-.00056 

3.00000 

.64048 

.20384 

.02819 

.14342 

-.00098 

3.10000 

.64229 

.19659 

.02859 

.14582 

-.00134 

3.20000 

.64397 

.18974  > 

.02896 

.14809 

-.00167 

3.30000 

.64554 

.18326 

.02930 

.15025 

-.00195 

3.40000  * 

.64700 

.17713 

.02962 

.15230 

-.00219 

3.50000 

.64836 

.17131 

.02991 

.15425 

-.00241 

3.60000 

.64963 

.16579 

.03019 

.15610 

-.00259 

3.70000 

.65082 

.16054 

.03044 

.15787 

-.00276 

3.80000 

.65194 

.15555 

.03068 

.15955 

-.00290 

3.90000 

.65299 

.15081 

.03090 

v  .16115 

-.00302 

4.00000 

.65398 

.14629 

.03111 

.16269 

-.00312 

4.10000 

.65491 

.14197 

.03130 

.16415 

-.00321 

4.20000 

.65579 

.13786 

.03149 

.16555 

-.00329 

4.30000 

.65662 

.13393 

.03166 

.16689 

-.00335 

4.40000 

.65740 

.13018 

.03182 

.16818 

-.00340 

4.50000 

.65814 

.12659 

.03197 

.16940 

-.00345 

4.60000 

.65884 

.12316 

.03212 

.17058 

-.00348 

4.70000 

.65951 

.11987 

.03225 

.17171 

-.00351 

4.80000 

.66014. 

.11671 

.03238 

.17280 

-.00353 

4.90000 

.66074 

.11369 

.03250 

.17384 

-.00355 
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APPENDIX  B 


PROGRAM  FOR  VALUE  OF  Srkj(a=L0,  oo) 
1.  PROGRAM(MATLAB) 


* 

fcSPHERE    (XLPHA-1.0) 

*  ! 

for  i-l:6; 

nu-i/10.;  i 

den**  1 5 .  *  ( 1-nu) ; 

sll(i)«(7-5*nu)/den; 

sl2(i)-(5*nu-l)/den; 

sl3(i)-sl2(i); 

s31(i)»sl2(i) ;  i 

s33(i)-sll(i) ; 

i 

%ALPHA«INPINITY  * 

% 

den2-8* (1-nu) ; 

ssll(i)-(5-4.*nu)/den2;' 

ssl2(i)--(l-4*nu)/den2;i 

ssl3(i)-nu/(2.*(l-nu) ) ; 

ss31(i)-0; 

ss33(i)-0; 

end ; 

resultl- [sll ; 812; sl3 ; s31 ; s33 ; ] ' 

result2»[ssll;ssl2;ssl3;ss31;ss33] ' 
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2.  Srkl  VALUE 


\         c 

Sll 

S12 

S13 

S31 

S33 

ru          ^s, 

re 

sultl   = 

C.l 

0.4815 

-0.0370 

-0.0370 

-0.0370 

0.4815 

0.2 

0.5000 

0 

0 

0 

0.5000 

0.3 

0.5238 

0.0476 

0.0476 

0.0476 

0.5238 

0.4 

0.5556 

0.1111 

0.1111 

0.1111 

0.5556 

C.3 

0.6000 

0.2000 

0.2000 

0.2000 

0.6000 

C.5 

0.6667 

0.3333 

0.3333 

0.3333 

0.6667 

re 

sult2    = 

0.1 

0.6389 

-0.0833 

0.0556 

0 

0 

0.2 

0.6563 

-0.0312 

0.1250 

0 

0 

C3 

0.6786 

0.0357 

0.2143 

0 

0 

0.4 

0.7083 

0.1250 

0.3333 

0 

0 

0.5 

0.7500 

0.2500 

0.5000 

0 

0 

0.6 

0.8126 

0.4375 

0.7500 

0 

0 
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APPENDIX  C 
PROGRAM  FOR  VALUE  OF  THERMAL  RESIDUAL  STRESS, 

cp  ct*  c-c  -i  ^  -l 

1.  PROGRAM(MATLAB) 

.r:q 

%This  program  computes  stress  values  and  constants 

X 

X  Initial  values 

nu-0 . 3 ; 

em«47.5e9; 

ei-427.0e9; 

sy-47.5e6; 

f-0.2; 

dt— 200 

ap-23.6e-6;  i 

as-4.3e-6; 

mu«0 . ; 

X 

X  Eshelby's  tensor  components 

sll-0.6786; 
s33»0.0000; 
S12-0.0357; 
sl3-0.2143; 
S31-0.0000; 

* 

X  Constants  f 

* 

app-(as-ap) *dt;  • 

g-em/(  (l.+nu) * (1.-2. *nu) )  ; 

h-ei/  { (l.+nu)  *  (1.-2. *nu) )  ; 

j-(em-ei)/((l.t-nu)Ml.-2.«nu)); 

a-f+sll*(l.-*); 
b-sl2*(l.-f) ; 
c-sl3Ml.-f); 
e-f+s33«(l.-f); 
d-s31*(l.-f); 
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continued 


tl-rl+r2+r3-ql+q2+q3 ; 

t2-(rl+r2-2*r3)/2.-(ql-q2+2*q3)/2.; 

t3«r4+r5+r6-p4-p5-p6; 

t4-(r4+r5-2*r6)/2.-(p4+p5-2*p6)/2.; 

t5»r7+r8+r9-pl-p2-p3; 

t6-(r7+r8-2*r9)/2.-(pl+p2-2*p3)/2.; 

xl-g* ( (l.-nu) «tl+nu* (t3+t5) ) 

x2-q« ( (l.-nu) •t2+nu« (t4+t6) ) 

x3-g* ( (l.-nu) *t3+nu* (tl+t5) ) 

x4«q*(  (l.-nu)  *t4+nu*  U2+t6) ) 

x5-q* ( (l.-nu) *t5+nu* (tl+t3) ) 

x6»q* ( (l.-nu) *t6+nu* (t2+t4) ) 

cl-( (x2+x4)/2.-x6)*(f/2.) ; 

c2-((xl+x3)/2.+x2+x4-x5+x6)*(f/2.)  ; 

% 

x  Temprature  dependent. parameter 

ep— c2*app/(2.*cl)+(l.-f)*sy/(2.«cl); 

% 

X  Elastic  moduli  of  the  matricies 

aaa-g* [ (1-nu)  nu  nu  0  0  0;nu  (1-nu)  nu  0  0  0;nu  nu  (1-nu)  0  0  0; 
0  0  0  2.*mu  0  0;0  0  0  0  2.*mu  0;0  0  0  0  0  2.*mu]#; 

* 

*        Plastic  strain  'in  the  matrix 

* 

eijp-[-ep/2  -ep/2  ep  0  0  0]' 

% 

\  Uniform  transformation  strain  of  inclusion 

eijts"[app+ep/2  app+ep/2  app-ep  0  0  0]' 

% 
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continued 


eijt-[ql* (app+ep/2)-q2* (app+ep/2)-q3* (app-ep)     . . 

p4* (app+ep/2)+p5* (app+ep/2)+p6* (app-ep)     . . 

pi*  (app+ep/2)+p2*  (app+ep/2)+p3*  (app-ep)    0   0   0] 
eij_c-[rl* (app+ep/2)+r2» (app+ep/2)+r3* (app-ep)     . . 

r4*  (app+ep/2)+r5*  (app+ep/2)+r6*  (app-ep)    . . 

r7* (app+ep/2)+r8* (app+ep/2)+r9* (app-ep)    0   0   0] 
% 

\  Average  internal  stress  in  inclusion 

% 

si j_i=aaa* (ei j_ c-eijt) 
% 

%        Average  stress  in  the  matrix 
% 
sijm—f  *sij_i 
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2.  CONSTANTS  OF  ELASTIC  MATRIX  AND  ELASTIC  INCLUSIONS 


B=Sn(l-f)  +  f 

D=SI3(l-f) 

F«S«(M) 

H=E1/(l+^)(l-2//) 

K=J((1-»/)B+i/C+i/F)-(1-i/)G 

M=J((1-i/)D+j/D+z/E)-*/G 

P=J((l-i/)B+i/C+i/F)-(l-i/)G 

R=J((l-^)F+^B+i/C)-^G 

Y=J((l-i/)B+2i/DHl-i/)G 

V=(QK-NM)/K 

X=(SK-RL)/K) 

Z=AH(1+//)(R-K)/K 

y'=(w-xV)/u 


c=s12(i-n  +  f 

E=S33(l-f)  +  f 
G=Em/(  1+i/)(1-2i/) 
J=(Effl-Ei)/(l+i/)(l-2J/) 
L=J((l-v)C+uB+uF)-vG 
N=J((l-4/)D+i/D+i/E)-^G 
Q=J((l-,/)D+//D+//E)-z/G 

S=J((l-i/)?+vC+vB)-vG 

U=(PK-NL)/K 
W=AH(l+i/)(N-K)/K 
Y=(TK-MR)/K 
X=(ZU-XW)/(UY-XV) 
z'^-YL-x'm-Ha'ci+^/K 
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APPENDIX  D 
FORMULA  OF  S1)ki  FOR  A  VARIETY  OF  INCLUSION 


1.  FORMULA  OF  Sikl(a=1.5  -  5) 


•Sill  -  S'lXU 


:\  «2  1 


8(1  -  Hit*--  I)      4(1  -*<,/ 


1  -  2*0  - 


9' 


4Ut*-  1) 


■ti 


Sxxu  - 


I 


2(1 -i/o)  I 

£  1 122  -  ^2211 


Cirt*-  l)       f  :»a2 

(I  -  2*0  +  — : — —  -  1 1  -  2*o  +  — — -k 


ia*-l)       I 


(«2-l)J 


1 


4(1  -*0)  12(«--  1) 


a*  3 

-(1-2*0) " A' 


»^'li:i3  -  ^'2233  -  ~ 


1 


tY* 


4(a--  1)    J 
1  f     .W 


2(  1  -  *0)  U»J  "  1)      4(1-  *0)  t(a-  -  1 ) 

-  (l-2*o) 


A' 


•*>'  !.'tll   ~  •*>.H22  s  ~ 


1 


2(1 -*0) 


1  -  2*o  + 


1 


la--  l)j 
1 


2(1  -*„)l 


3 

1  -  2*n  + ; \t 

2(«--  i»r 


where  *o  is  Poisson's  raiiu  of  a  matrix,  c»  is  aspect  ratio  ol'  a  fiber 
(  =  //d),  and  g  is  given  by 


5 


nr 


-Ima2-  l)w-cuih-'  «J 
(a* -I)*2  r 
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2.  FORMULA  OF  S"kl(a=l) 


q  7  -  bv 

Sn..=  15(l   -i>) 


q  5y  -  1 


c  5^  —  1 

b..33"15(i    "    I/) 


c  ^  -  1 

b22l."15(l     "    V) 


o       _      7  —  bv 

52222"  15(1    -  V) 


c       _    5^  —  1 

2233~15d     -    v) 


o  5^  —  1 

S33.rl5(l    -    I/) 


c        _       bv  —   1 
3322~15(1    -    I') 


S  7  -  5^ 

3333       15(1    ~  ") 
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3.  FORMULA  OF  S-yC^oo) 


c  5  -  4v 

Su..-8(1  =T) 


Q  k   -    1 

1122~  8(1    -    V) 


Sim"^-  ") 


S22..=  8(l    -    U) 


S       _    5  -  4t/ 

2222       8(  1    —   Iy) 


2233"  %  1    "    *) 


s     = 

3311 

0 

s     = 

3322 

0 

S         = 
3333 

0 
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APPENDIX  E 
PROGRAM  FOR  CREEP  DEFORMATION 


PROGRAM(MATLAB) 


ep0=0.0Gl; 

beta=C . 0465; 

k=7.2e-7; 

for    tt=l:91; 

kv(tt)=tt; 

c-1.0«10*( { (tt-l)/30.) )  ; 

eep(tc)=epO«  (1  . +beta*t~  (1/3)  )  »exp(k*t)  ; 

end; 

plot(kv,eep) 

seir.ilogx; 

xlabel ( 'Time    t '), ylabel (' Strain   eps * ) 
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